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Abstract 

We introduce a method for the coherent acceleration of atoms trapped 
in an optical lattice, using the well-known model for stimulated Raman 
adiabatic passage (STIRAP). Specifically, we show that small harmonic 
modulations of the optical lattice amplitude, with frequencies tuned to the 
eigenvalue spacings of three "unperturbed" eigenstates, reveals a three- 
state STIRAP subsystem. We use this model to realize an experimentally 
achievable method for transferring trapped atoms from stationary to mo- 
tional eigenstates. 

1 Introduction 

Stimulated Raman adiabatic passage (STIRAP) is a method for achieving co- 
herent transitions between quantum states by applying two coupling fields in a 
non-intuitive pulse sequence. The frequencies of these fields arc tuned to the 
eigenvalue spacings between the "initial" and "target" states and a third "in- 
termediate" state. When the coupling of the target and intermediate states 
precedes that of the initial and intermediate states in time, population transfer 
from the initial to target state is achieved. In the adiabatic limit, the transition 
is 100% efficient and involves no occupation of the intermediate state. The use 
of STIRAP for atomic and molecular systems was first demonstrated experi- 
mentally by Gaubatz and coworkers [H, [3| , who achieved population transfer 
between vibrational levels in a beam of sodium molecules. Further references 
on the STIRAP transitions in atoms and molecules can be found in the review 
by Vitanov et al [3]. In atom optics experiments, STIRAP has been used for 
coherent momentum tranfer 0, Q , and velocity-selective coherent population 
trapping for laser cooling of trapped atoms [g, J\ . Extensions of STIRAP, with 
a particular focus on the influence of quantum chaos, have been studied by Na 
and Reichl fs','^. Na et al [l3| have also use STIRAP to control the isomerization 
transition of HOCl. 
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Figure 1: The 3-level "ladder" STIRAP system (Eq. [1]) with parameters 
A = 0.05, Wi — W2 = 1 and h — 2. Couphng fields are applied with fre- 
quencies equally detuned from the spacings of the unperturbed energy levels 
by A (a). Adiabatic variation of amplitudes of the coupling fields (b) in the 
manner described by Eq. ^ affects a transition of the |1) eigenvector between 
basis states \a) and |6) (d). The eigenvalue corresponding to this state remains 
unchanged at zero throughout the transition (c) 



As a theoretical model, STIRAP can be defined by the adiabatic behavior 
of a three-level system, which in some basis (|a), |c)) is represented by the 
Hamiltonian 

. / Wiit') \ 
H{t')^-- W,{t') -2 A W2{t') , (1) 
V ^2(i') / 

under the variation of the parameters Wi and W2, with A a constant. This 
model, which we refer to as the "STIRAP model" hereafter, was first introduced 



by Kuklinski et al , following significant work by Hioe et al [ij, ll3|, , to 
succinctly describe the experimental results of Gaubatz ei a/ [ij . Starting from 
a system in which the pairs of states (|a), and (|6),|c)) are each dipole- 
coupled by monochromatic electric fields, the STIRAP model is derived by 
applying the rotating-wave approximation and assuming an equal detuning of 
the coupling frequencies A (see Figure [1]). In that particular system the Wi{t') 
are the Rabi oscillation frequencies corresponding to the two couplings. The 
STIRAP transition, however, does not depend on the physical system from 



which Eq. ([T]) is derived. In a novel appHcation of this model by Eckert et al 



15|, for example, the functions Wi were related to the spatial separations of 
three optical microtraps in order to induce coherent transport of atoms between 
the ground states of the two farthest separated traps. 

The matrix in Eq. ^ allows for a transition of the type described above 
because of the existence of the eigenvector 

\l{t')) = cose{t')\a) ~ sme{t')\c) , (2) 
where ta.nd{t') = . Under the conditions 

lim }. — + and lim } ^ , (3) 

the adiabatic evolution of state |1) is from \a) to |c). Thus, the transition is 
achieved by first coupling the upper two levels and then coupling the lower two 
(in some continuous, e.g. Gaussian, manner). Moreover, because of the form of 
the eigenstate the state \b) remains unoccupied throughout the transition. 

In this paper, we apply the method of STIRAP to the motional states of 
atoms in an optical lattice. As was first shown by Graham, Schlautmann and 
Zoller the interaction of a single transition in an alkali atom with a pair 
of counter-propagating lasers can be reduced to an effective Hamiltonian for 
the center-of-mass motion of the atom in a cosine potential. Modulation of the 
laser amplitudes and/or the introduction of laser pairs with offset frequencies 
introduces a periodic time-dependence (see Appendix A). Here we will analyze 
the "two-resonance" system described by the effective Hamiltonian 

Ho{t) — + i^Q [cos X -f cos(i; — ujat)] , (4) 

where each of the two cosine terms are produced by a pair of counter-propagating 
lasers with kq proportional to the square of the laser amplitudes. We show, 
using perturbation analysis of an associated Floquet Hamiltonian, that small 
harmonic modulations of the laser amplitudes can be used to affect a STIRAP- 
like transition from a state localized in the stationary cosine well into a state 
localized in the traveling cosine well. 

In Section [2] we present an analysis of the time- independent "quantum pen- 
dulum" system which reveals a three-level STIRAP model in the regime of small 
perturbation. Section [3] contains a slightly modified method in order to obtain 
a STIRAP-like model for the time-dependent "two-resonance" system. In each 
case, adiabatic results are compared to numerical evolution of the Schrodinger 
equation and transitions of nearly 100% effeciency are observed. Concluding 
remarks are presented in Section U) 



2 STIRAP transitions in the Quantum Pendu- 
lum 

Before analyzing the two-resonance effective Hamiltonian presented in the in- 
troduction, we will consider STIRAP-like transitions within the quantum pen- 
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dulum system 

^pend 



Hpend ^ + Ko COS X . (5) 



This is the simplest type of effective Hamihonian for optical lattice experiments, 
achieved with a single pair of counter-propagating lasers with equal frequencies. 
The parameter kq is proportional to the square of the electric field amplitude. 
Because of experimental techniques which can limit momentum values to the 
integers (see Appendix[X| , the eigenstates can be considered considered spatially 
periodic with period 27r. The position space solutions to the eigenvalue equation 
HpendlXn) = En\xn) are the Mathieu functions {x\xn) (n S Z) with n even 
labeling even-parity functions and n odd labeling odd-parity functions. 

To affect a transition in this system using STIRAP, we add a time-periodic 
modulation of the lattice amplitude of the form 

XV {t) = A cos X [ki cos(fiit) + K2 cos(r22i)] , (6) 

where A is small, and fli and are commensurate with §^ = ^ (jrti G Z). 
It is useful to write Qi and Q2 in terms of a common frequency lu such that 
fli = mi uj and = m2 ui, where u = ^ and T is the periodicity of the 
perturbation. This perturbation is achieved experimentally by modulating the 
intensity of the counter-propagating laser radiation about the kq value, meaning 
that the lasers' electric field amplitude E{t) should take the form 

^ Ko + A [ki cos{miUJt) + K2 cos{m2UJt)] . (7) 

For now, we will consider the coefficients ki and K2 to have constant values. 
Perturbation analysis will reveal that, for small values of A, there exists a 3-state 
subsystem identical to the STIRAP model, parameterized by these coefficients. 
A STIRAP-type transition can then be affected by adiabatic variation of ki 
and K2 in the manner described in Eq. The justification of this time- 

parameterization of the Ki, following a Floquet analysis of the system where 
they are considered constant, is provided in Appendix [Bl 

Having required that the frequencies fii and Q2 are commensurate, we can 
analyze the dynamics of the full system 

= [Hpend + XVit)] (8) 

using Floquet theory, which we review briefly here. Any solution of Eq. ((8]) can 
be written in the form 

IV'.W) =e^"°*l0aW>, (9) 

where the Floquet eigenstate \(f)a(t)) is periodic in time with period T, and ea 
is called the Floquet eigenvalue. Plugging this solution into the Schrodinger 
equation we arrive at the eigenvalue equation 



d 

Hpend + XV{t)-i — 



(10) 
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where the Floquet Hamiltonian Hp is a Hermitian operator in an extended 
Hilbert space which has time as a periodic coordinate [3, Hj 01 • Diagonaliza- 
tion of Hp in some appropriate basis in this space yields the Floquet eigenstates 
and eigenvalues. The Floquet eigenstates have an infinite multiplicity with re- 
spect to the solutions of the Schrodinger equation in the sense that for any Flo- 
quet eigenstate |0) with eigenvalue e, there exists another solution |0) = e"^"*|(/)) 
{Q e Z) with eigenvalue e = e + Qoj. These two Floquet states are, however, 
associated to the same physical state, i.e. = e~"^*|0) = e"**^*]^). The impli- 
cation of this fact is that the dynamics of the physical system can be understood 
by considering only those Floquet states whose eigenvalues appear in a single 
"zone" e* < e < e* -I- w, labeled by the constant e*. 

The STIRAP model system is derived by applying perturbation theory to Eq. 
ifTO]) . where the two frequencies in V{t) are chosen to "couple" three pendulum 
eigenvalues, at a particular value of kq, in the manner shown in Figure [T^. 
Although, in general, the ratio of these eigenvalue spacings is not rational, the 
equal detuning A allows for fli and to be chosen as commensurate. More 
precisely, given three eigenvalues of the quantum pendulum Ea < Ei, < Ec, any 
pair of integers (mi, TO2) uniquely determines A and w via the coupled equations 

miuj = Eb - Ea - A ^^^^ 

7712 LO ^ Ec — Ef, + A . 

Eliminating uj we obtain an expression for A 

^^m2^i-mi^2^ (12) 
mi + m2 

where we have defined oji = E^ — Ea and L02 = E^. — Ei,. We can see that the 
integer vectors m = (1711,1112)^ which minimize A are those which satisfy 

m • i7 w with 1^ = {—LU2,i-ui)'^ (13) 

Therefore, the best choices for m are those for which mi/m2 are the best rational 
approximants of the ratio w = U!i/uj2. 

Let us denote the unperturbed Floquet Hamiltonian {Hp when A = 0) as 

H"p = Hpend - . (14) 

In the extended Hilbert space. Hp has normalized eigenvectors of the form 

{t\Xn,q)^{t\q)\Xn)=^e^''''\Xn) , Q^Z. (15) 
V 1 



The corresponding eigenvalues are e° g = En + quj. We are interested in the 
dynamics of an initial population of atoms localized in the pendulum state \xa)- 
With A = 0, state \xa) is represented in a particular zone of Floquet eigenvalues 
of Hp by a Floquet eigenstate |xa,9a) with eigenvalue ea.q^ in that zone. The 
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Figure 2: The lowest seven eigenvalues of the quantum pendulum (a) as a func- 
tion of the parameter kq (solid lines are the eigenvalues of even-parity states, 
dashed odd). The coupling of levels Eq-E^ and E4-EQ at = 8 via the per- 
turbation V{t) requires a Floquet treatment of the pendulum system. The 
corresponding Floquet eigenvalues (for oj = 0.5235) are shown in (b). The de- 
generacy and near-degeneracy of the the three coupled eigenvalues can be seen 
in the enlarged (inset) view. 
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coupling frequencies Qi and have been chosen in Eq. ([TT|) such that e° , the 
Floquet eigenvalue in that zone corresponding to the physical state |xc), is equal 
to eg and the eigenvalue eg is offset from this value by A. The degeneracy 
of these two eigenvalues and the near-degeneracy of the third requires that any 
perturbation analysis must be performed in the degenerate form [20|. Therefore, 
we expand the extended Hilbert space state {(j)) in powers of the small parameter 
A 

\cj)) = |(/.(o)) + A|0(i))+A2|0(2)) + ... (16) 

and take the zeroth-order state to be a superposition of the three degenerate or 
near-degenerate eigenstates of Hp : 

10^"^) = CalXa, qa) + Cb\Xb, Qb) + C,\Xc, Qc) ■ (17) 

The eigenvalue is likewise expanded in orders of the small parameter 

e = e(") + Ae(i) + A^e^^) + . . . . (18) 

For brevity of notation, we write the unperturbed eigenstates \a) = |xa,9a), 
\b) = \Xb, qb) and |c) = |xc, Qc), with associated eigenvalues eg = eg eg = eg 
and eg = eg ,^^, respectively. 

Retaining terms up to first-order in A, we obtain an isolated subsystem of 
Equation (fTT ^ 



(19) 



where e' = eg -I- Ae*^^-* and we have assumed that A is of order A. The matrix 
elements of the perturbation are determined in the following way: 

/7r /u 
{(it\t) {t\q3){xt\cosx\xj) 
-7r/u 

X [ki cos{miUjt) + K2 cos{m2U)t)] dt 

_ {xi\cosx\xj) 



eg + XVa.a 


XVa,b 


XVa.c \ 


/ Ca 


XVb.a 


eg + XVb,b 


XVb,c 


Cb 


XVc,a 


\Vc,b 


eg + AK,c / 






2 



(20) 



The selection of fii and il2 guarantees that the q indices for these states satisfy 

qb - qa = mi and qc - qb = m2 ■ (21) 

Therefore, the only non-zero matrix elements of the perturbation are Va^b = 
and Vb^c = Vc^b- Further, we can write 

Va,b = Vb^a = Kl Va,b and Vb^c = Vc,b = K2 Vb,c (22) 
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Figure 3: The Husimi functions of the first four even-parity eigenstates of the 
quantum pendulum on classical phase space (xo,po)- The states are, clockwise 
from top left: |xo),IX2),IX6),IX4)- 
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where the numerical values of 

_ {X^\cOSx\Xj) , . 

are calculated using the Mathieu functions. Subtracting e° {Ca,Cb,Cc)^ from 
both sides of Eq. (flQ)) and redefining e' — ^ e', we arrive at 

>^KlVa,b 

A Ki Vb^a A A K2 Vb,c | | Cfc | = c' | Cfc | , (24) 

^K2Vb.c 

which is equivalent to the STIRAP model system. Thus, in the limit of small A 
(and A), the parameters ki and K2 can be adiabatically varied in the manner 
described in the introduction to affect a transition between the unperturbed 
states I a) and |c). 

We now provide a concrete example on which to demonstrate the analysis. 
Figure [2^ shows the energies of the first few eigenstates of Hpend as a function 
of kq. Husimi representations 13, 21| of the even eigenstates, at kq = 8, are 
shown in Figure [3l At this value of kq, we choose energies Ea = Eq, Eb — E4 
and Ec — Eq for coupling. These energy levels have spacings u>i — Eb — Ea = 
12.5668395 and uj2 ^ E^ - Eb ^ 3.6630472 with ratio 

w = 3.43070639... = [3, 2, 3, 9, ...] = 3 + — . (25) 

^ + 3 I 1 

Thus, the best rational approximates of w, found by truncating the continued 
fraction, are {j, |, ...}. Using the third approximation, the modulation 

frequencies shown in the example have been chosen to be SI 1 = 24a; and ~ 7uj, 
giving uj = 0.5235 and A = 1.766 x 10"^ 

In Figurej^bi the Floquet eigenvalues of Hp are shown in the zone e* = 0. In 
the inset figure, one can see that = £9,12 and e° = Cg^-ig are equal at acq = 8 
and e" = £4-^12 is offset by A. Using the Mathieu functions, the perturbation 
matrix elements are calculated to be 

Va^b = Vb,a = Kl Va,b = -1-16 X IQ-^^^i (26) 

Vb,c = Vc,b = K2 Vb^c = 2.50 X 10-1^2 . (27) 

To accomplish a STIRAP transition from \a) to |c), we give ki and K2 
Gaussian dependence on an adiabatic time parameter t' (see Figure H^): 



Ki{t') — exp 



(<' - U) 



21 



2a 



(28) 



The conditions of Eq. ^ are satisfied by setting ti = —t2 = 1.0 and ai = (T2 = 
1.0. Figure |4] shows good agreement between the adiabatic dynamics of the 
model system in Eq. (|24)) (|4)d and d), and that of the full Floquet Hamiltonian 
([it and e). 



9 




Figure 4: Adiabatic evolution of the pendulum system under STIRAP coupling, 
comparing the dynamics of the model in Eq. (b and d) to that of the full 
Floquet Hamiltonian in Eq. ([8]) (c and e) for A — 0.1. The coupling field K2 is 
seen to have a greater effect on the eigenvalues at t' = — 1 than that of ki at 
t' — I, since Vb^c » Va,b- This asymmetry also shifts the transition of the |1) 
eigenstate from pendulum states \a) to |c) to a later time (d and c). 
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The implementation of this transition in an experimental system (or the 
numerical evolution of the Schrodinger equation) is not dependent on the time- 
periodicity which we have required thus far. Floquet analysis has proven an 
essential theoretical tool for revealing the existence of the STIRAP model, but 
the method has introduced no upper limit on the integers mi and m2 whose 
ratio approximates LOxjiOi. Therefore we may choose the coupling frequencies to 
be resonant (fJi = iO\ and ^2 =1^2) to any desired accuracy. The results for the 
numerical evolution of the effective Schrodinger equation ([5]) are shown in Figure 
[5^, for the case of both resonant and near-resonant coupling (A — 1.77 x 10^"^) 
with A = 0.1. The evolution was performed over a set time period [0, itot] 
with initial condition |(Xa|V')P = li Gaussian parameters for the coupling 
amplitudes ki and K2 of ui = (T2 = O.lttot, ti = 0-6ttot and t2 ~ OAttot- It 
is seen that resonant coupling provides a more rapid approach to the adiabatic 
behavior. In Figure [Sb, good agreement is seen between the resonant evolution 
of the full effective Schrodinger equation and the adiabatic predictions of Figure 

m 



3 STIRAP transitions from stationary to mov- 
ing atoms 

We now consider the case in which the "unperturbed" Hamiltonian has the form 
of Eq. (HI, which consists classically of a stationary cosine wave and a cosine 
wave that travels through phase space with a speed ujq. Our goal is to cause a 
coherent transition of an entire cloud of trapped atoms from a state localized 
in the stationary wave (in the sense of its Husimi distribution) into a state 
localized in the traveling wave, so that the entire collection of atoms changes 
velocity from v — to v — ujq. 

Our approach is analogous to that of the previous section. We apply per- 
turbation theory to a Floquet eigensystem of the form 

HF{t)m)) = [HUt) + >^v{t)\ m)) = em)) , (29) 

where Hp is periodic in time with period T = 2t:/uj and the perturbation 
operator V has the same form as in Eq. ([6]). Again we find that, in the limit of 
small A, there exists an isolated three-level subsystem in which a STIRAP-type 
transition between eigenstates of Hp can be induced. The construction of Hp, 
however, differs significantly from the previous section because of the explicit 
time-dependence of the two-resonance Hamiltonian. In the pendulum analysis 
the frequencies fli and fl2 were chosen to couple the energies of pendulum 
eigenstates. Here, these frequencies are chosen to couple the eigenvalues of the 
two-resonance Floquet Hamiltonian, 

d 

Hpit) =p'^ + Ko [coax + coa{x - ujot)] ~ ; (30) 

within a particular zone. Selection of coupling frequencies such that they and 
0^0 are commensurate allows for Floquet analysis of the full system, but requires 
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Figure 5: The numerical evolution of a state initially localized in the \xa) pen- 
dulum state, under STIRAP coupling of the energies. The probability of a 
non-adiabatic transition (transition to any state other than |c)) is plotted ver- 
sus the width of the gaussian field a ^ ai = (a), for resonant (triangles) and 
near-resonant coupling (squares) with A = 0.1. The solid lines are the values 
predicted by evolution of the three-state model (see Appendix [B] for a justifica- 
tion of this use of the model). Figure (b) shows the numerical evolution, under 
the effective Schrodinger equation, of a state classically localized within the pen- 
dulum well to one localized on the separatrix (axes on inset Husimi functions 
are the same as in Fig. [3]). This evolution corresponds to the point at cr = 6000 
with resonant (A = 0) coupling in (a). 
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that the eigenvectors of H'p be translated from their natural Hilbert space, 
containing functions periodic in time with period T = 27r/a;oi into the space 
containing T-periodic functions of time. The relevant Floquet eigenvalues asso- 
ciated to the eigenvectors in this latter space take near-degenerate values and 
perturbation analysis leads to similar results as the previous section. [It should 
also be noted that although we call the two-resonance system in Eq. (|5D|) an 
"unperturbed" Hamiltonian, it is not analytically solvable. Perturbation the- 
ory will be a useful tool to demonstrate the existence of a STIRAP-like model 
for this system, but the eigenvectors of Hp and all related quantities (e.g. the 
matrix elements of the three-level system) must be determined numerically.] 

We begin by constructing Floquet eigenvectors |^°) of Hp (the "overbar" 
will be used to indicate that these vectors belong to the Hilbert space 7Y, defined 
below). We will assume that the parameters kq and have constant values, 
which may be set arbitrarily. The only limitation on this choice is that, given 
070, Ko should be chosen such that the set of eigenvectors with eigenvalues in 
a particular zone contains one state localized purely in the stationary cosine 
wave and one state localized in the traveling wave (i.e. a kq value far from 
avoided crossings involving the eigenvalues of these states). The eigenvectors 
|^°) lie in the extended Hilbert space 7i = 8 (g) T, where O is the space of all 
27r-periodic, square-normalizable position-space functions and T is the space of 
all T-periodic, square-normalizable functions of time. We select the complete 
set of momentum eigenstates \n) (see Appendix A) as a basis in G and the 
analogous eigenstates \q) as a basis in T, yielding normalized basis vectors in 
the extended space which can be written 

{x, q) = {x\n) {t\q) = -^e^^e^^""* , (31) 
with n^q eIj. The eigenstates \<j)^) can then be written 

\~^°{t)) = E ^e^''"°*(n,g|4°) \n) , (32) 

n,q ^ ^ 

where \4'a{^)) — {A4>a) the coefficients (n, are determined by diago- 

nalization of Hp in 7Y. 

We select a zone e* < < e"^ -I- ojq within which to perform a coupling of the 
eigenvalues e° of H^. Two of these eigenvalues, denoted e° and e°, are those of 
the states localized in the stationary and traveling waves, respectively. A third 
eigenvalue is chosen with the restriction that the corresponding eigenvector 
is localized "nearby" in phase space (the matrix element of cos x between this 
and the other two vectors should not be vanishingly small). As before, the 
coupling frequencies Sli and ^2 must be chosen to be commensurate. In this 
case, however, analogous equations to Eqs. (fTTj) cannot be solved simultaneously 
with the requirement that wq is likewise commensurate: 

'jjQ = mQLo (too G Z). (33) 
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Therefore, in the following, we will relax the constant detuning requirement and 
allow for two independent detunings defined by the equations 

Given any integer vector m = (toq, mi, m2)'^ , Eqs. ([55)1 and ([M)) can be solved 
for (Ai, A2,w). Eliminating uj and defining uji = l^ — and uj2 = — e^, we 
obtain 

(35) 

A2 — . 

mo 

Thus, we see that the integer vectors m which simultaneously minimize the two 
detunings will be those closest to the vector perpendicular to the plane defined 
by = (wi,— wqiO)'^ and iJ^^^ = (— W2, 0, wq)"'"- This perpendicular vector 
is of course ft = (wq, ti-"!, 0^2)'^, and the problem of minimizing the detunings is 
reduced to finding the best integer approximate of n or, equivalently, finding 
the simultaneous pair of rational approximants for (wi/wq, W2/a;o)- 

In the context of the three-level model system presented in the introduction, 
non-equal detuning of the coupling frequencies leads to a Hamiltonian of the 
form 

/ Vt^i \ 

H = -- \ Wi -2Ai W2 . (36) 

2 V W2 A2 - Ai / 

It was recognized by Kuklinski et al [11] that this system could allow for a 
STIRAP transition, despite the absence of an analytical result analogous to Eqs. 
(UJ and (lH), as long as the condition ^/Wf+Wj » IA2 - Ai| is satisfied. In 
order to satisfy this requirement, and that of small perturbations, we will seek 
integer vectors rh which provide detunings IA2 — Ai| << Ai << 1. 

Since the full, perturbed system Hpit) is periodic in time with period 
T = 27r/a;, we must determine the eigenstates of the unperturbed Floquet 
Hamiltonian in the extended Hilbert space TC = Q ^T, where T is the Hilbert 
space of T-periodic functions. These eigenstates can be expanded in H as 

n,q ^ 

where the g-eigenvectors now have the time-periodicity of T. Since the Schrodinger 
equation for the unperturbed system 

d 

^^IV^i)) = {p^ + Ko[cosx + cos{x~ujot)]} Hit}) (38) 

can be viewed as time-periodic with either period T or T = moT, a physical 
solution \ilja{t)) can be written, using Eq. ((9]), in terms of a Floquet state with 
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either periodicity. Equating these two representations, we obtain a relationship 
between the Floquet eigenstates of in spaces and Ti: 

exp [-^llt] \4>lit)) = A exp [-^elt] |<^^(t)) , (39) 

where e° is the eigenvalue associated to 10° (i)) and A is a proportionality con- 
stant. Equating coefficients of the momentum eigenstate \n) in Eqs. (I32p and 
((571) . we find 

exp[z<7moc^i] (n, = J] A' exp [z [q'u; - + 6°) (n, q'\<Pl) , (40) 

where A' is again a constant. A non-trivial solution to this equation requires 
that the eigenvalues satisfy e° — — Quj, where Q is an integer. We see, then, 
that associated to each Floquet eigenstate in 7i is a family of eigenstates in H: 

{\€,q).^Iq} QeZ, (41) 

with Q ~ ^ Q^- Selecting a particular value of Q, Equation becomes 
exp [iqmaut] {n, q\^°) = ^ A' exp [i{q' ~ Q)iut] (n, <z'|<^^i,Q) . (42) 

q q' 

Equating coefficients of the exponentials, we find 

whence Z 

{n,q\€,Q)^{ ^ (43) 

otherwise , 

where we have set A' ^ 1 under normalization. Therefore we see that the 
unperturbed eigenstates in the space H have non-zero coefficients (n, g) 
only at mp-separated values of q, with an offset of Q from q — 0. 

Within a particular zone, we denote the unbarred eigenvalues corresponding 

to {eO,eO,e°} as [4,q^^4.q,^^IqJ\^ with values related by 

and corresponding eigenstates 1 1 ) : I Qj, ) , I q,, ) | • The Q-indices of these 
states are related by Qa~ Qb = mi and Qb — Qc — 

Perturbation analysis of Eq. (|29p is now performed by expanding the eigen- 
state and eigenvalue e in powers of A. Assuming that \4'aQ^) is initially 
occupied with probability one, and taking into account the near-degeneracies of 
Eq. p4)) , the zeroth-order term in the expansion of the perturbed eigenstate is 
chosen to be of the form 

1'/'^°^) = CalCgJ + ^^^>I<qJ + CM,qJ ■ (45) 



15 



Retaining terms up to first order in A and making the assumption that Ai and 
A2 are of order A, we obtain 

e° + Wa.a AK,b \Va,c \ ( Ca\ f Ca\ 

XVb.a eg + AH,6 XVb,c Cb = e' C'b , (46) 

XVc.a XVc.b e2 + AK,c J \Cc J \Cc I 

where e' = eJJ + Ae'-'^' and the matrix elements are calculated, defining |i) = 
\4>i,Qi), as follows 

= ((zlfb')) = mn,q){n\ cos x\n'){n',q'\3)] (47) 

n.q^n' .q' 

X \^ni{q\cos{miLLii)\q') + H2{q\cos{m20ji)\q')\ 
= i ^ Ki(n|cos£|n') [(^|'^,9)(?^^g + ml|j) + (i|n,g)(7i'g-mi|j)] 

n,q,n' 

+ K2(ri| cosi|7i') [{i\n,q){n' ,q + m2\j) + {i\n,q){n' ,q - m2\j)] ■ 

Recalling the structure of the states \<i>iQ.) given in Eq. (|43l) . we see that the 
sum 

^(i |7i, q) {n',q + m\j) = Y,{4'Iq, \n, q) {n',q + m|0O ) , (48) 

can be non-zero only when Qj+m — Qi + k mo with k Cz Z. Thus, the only 
non-zero matrix elements of V in Eq. I|46p are 



Va.b = ^ ((a|n,(7)(n + 1, g - mi|6) + (a|n, - 1, 9 - r7ii|6)) 



11 
4 

n,q 
n.q 

^ ((4°|(cos^ 01)10°)) ^^g^ 

and 

Vb,c = -^X! ((^ki'9)("- + ^,q- ™2|c) + {h\n,q){n -l,q- m2\c)) 

n,q 

= ^J2i(^S\n,q){n+l,qm + {^S\n,q){n~l,q\^^)) 

n,q 

((0°|(cos^ 01)10,")) 
= '^2 ^ (50) 

where the second and third equalities for each matrix element have been written 
in terms of the Floquet states in H, using Eq. (|43)l . Subtracting e°(Ca, Cfc, Cc)^ 
from both sides, redefining e' — e° — > e', and defining Va^b = V^j^b/Ki and Vb.c = 
yb.c/1^2, Eq. becomes 

Aki Va,b 

Aki Vb,a Ai Ak2 U6,c II Cb I = e' I Cb | , (51) 
Ak2 Vnb Ai - A2 
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which is equivalent to the desired model system in Eq. 

Again, we provide a example system on which to demonstrate the analysis. 
Figure [5^ shows some eigenvalues e° of the two-resonance Floquet Hamiltonian, 
in the zone labeled by e* = —u}q/2 with = 6.180339887, plotted as functions 
of Ko- A triplet of eigenvalues e° < e° < e° has been chosen at kq = 1 for 
STIRAP coupling. The Husimi functions of the three corresponding eigenstates 
are shown in Figure [T] The values of these eigenvalues satisfy loi = e ° — e ° = 
1.67227495 and uj2 = e° - = 1.14207065. Therefore we seek simuhaneous 
rational approximates (toi/too, m2/mo) to the pair (0.270579771, 0.184790913). 
Performing a numerical exhaustive search, we find that the integer vector to = 
(325,88,60) provides detunings Ai = -1.17 x lO'^ and A2 = -1.08 x lO^^, 
which satisfy the required conditions of Ai ~ 0[\) « 1 and |Ai — A2I << A. 
The unbarred eigenvalues in the zone e* = —uj/2 — —ujq/{2 x 325) are shown 
in Figure [6b. The detunings Ai and A2 can be seen in the enlarged section 
of the graph (inset), separating e° _gg from e° and 6° -1291 respectively. The 
coefficients 

0.4442 

Va,b = (52) 
-0-0673 

Vb,c = ^ (53) 

are calculated using Eqs. and ([50]) after numerical diagonalization of Hp 
in H. 

The adiabatic dynamics of the model system in Eq. ([ST]) and the full system 
in Eq. ()29|) are shown in Figure [8l using the same parameterization of the Ki 
as in Eq. ([28]) and A — 0.02. Although good agreement is seen between the 
two, it is evident that the STIRAP transition between eigenstates \a) and |c) 
is not achieved in either case. The reason for this failure is a narrow avoided 
crossing at t' « —2.5 between the eigenvalues of adiabatic states |1) and |3) (see 
inset of Figure [5^), which affects a transition between unperturbed states \a) 
and |c) before the STIRAP transition. This type of avoided crossing, reversing 
the effects of the desired transition, will always exist in the adiabatic limit when 
a matrix of the type given in Eq. is used for STIRAP evolution because of 
the non-degeneracy of the eigenvalues of \a) and |c). Although this model does 
allow for a broad STIRAP-type transition, the resulting change in character of 
the adiabatic state |1) as t' passes from —00 to 00 requires that its eigenvalue 
change from to Ai — A2 . 

The problem with the adiabatic model can be avoided in the numerical or 
experimental achievement of a STIRAP transition in one of two ways. First, it is 
possible to achieve a non-adiabatic evolution of the system which is slow enough 
to guarantee a STIRAP transition, but too rapid to "see" the problematic sharp 
avoided crossing. Second, one can abandon time-periodicity and apply resonant 
coupling fields, reducing the model to the classic form in Eq. ([T]). The efficacy 
of both methods can be seen in Figure [9^, where long-time evolution of both the 
model (see Appendix[B| and the full Schrodinger equation yields a STIRAP-like 
transition in the case of resonant and detuned coupling fields. As in the pen- 
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Figure 6: Eleven eigenvalues of the two-resonance Floquet Hamiltonian with 
loq ~ 6.18, viewed as an operator in 7i, plotted as a function of kq (a). The 
corresponding "unbarred" Floquet eigenvalues of the same Hamiltonian, viewed 
as an operator in the space Ti. with lo = a;o/325, are shown in the zone labeled by 
e* = — ti>/2 (b). The near-degeneracy of three eigenvalues , e° and e° g^ 
at Kq = 1 can be seen in the inset. It is evident that there are other eigenvalues 
nearly degenerate with these three, however these need not be considered in Eq. 
(|45|) since their respective Q-values will yield zero-valued matrix elements. 
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Q. 

1 ■ 



-1 

_2l ^ , ^ 

2 4 6 

Figure 7: The Husimi representations of three Floquet states of the two- 
resonance Hamiltonian {\<f>a{t)), \4'b{t)) a-nd clockwise from top left) 
viewed at time t — 0. The parameter values are kq = 1 and loq ~ 6.18. 
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Figure 8: The adiabatic dynamics of the 3-level model in Eq. (|5T|) (a and c) 
and the three relevant states of the corresponding full two-resonance Floquet 
system with ^ 6.18 (b and d) for the example at kq = 1. The perturbation 
amplitude functions Ki{t') and K2{t') are the same as in FigureH^, and A = 0.02; 
model parameters are as determined in the text. The eigenvalues under the 
application of the perturbation are show in (a) & (b). The influence of Ki(t') on 
the adiabatic eigenvalues is stronger than that of K2{t') because \va,b\ > l^b.cl- 
Overlaps of the |1) adiabatic eigenvector with the unperturbed states are shown 
in (c) and (d). The inset in (a) shows a sharp avoided crossing which prohibits 
the STIRAP-like transition in the adiabatic limit. 
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dulum case, the resonant coupling provides a faster approach to the transition. 
The evolution of a state initially prepared in the stationary wave eigenstate \a) 
is shown to pass into the traveling wave state |c) under resonant coupling in 
Figure [5)3. 

Using the non-dimensionalization presented at the end of Appendix \^ we 
can return to dimensional variables and determine the experimental conditions 
necessary for a STIRAP transition of the type described here. Consider a system 
of cold cesium atoms interacting with a system of lasers tuned near the D2 
transition (as in References [12, [Ij]), yielding a recoil frequency of ujr ~ 1.3 x 10* 
Hz. For the example considered above, the traveling cosine wave is therefore 
generated by counterpropagating lasers with frequencies offset by 6uj/2tt « 10 
kHz; amplitude modulation frequencies corresponding to f2i/27r and r^2/27r are 
11 and 7.7 kHz, respectively. The relationship between dimensionless time t and 
physical time tphys for this system is 

tphys ^ l-O ^ 10-^ t . (54) 

Therefore, the near-100% transition between the stationary and travelling lat- 
tices shown in Figure ^jp would require half a second in the laboratory. The 
approach to adiabatic evolution can be achieved more rapidly by increasing 
the coupling strength A [l3|, as long as the assumption of small perturbation 
(A << kq) remains valid. In numerical experiments, we were able to decrease 
the transfer time for the preceding example (with resonant coupling) by a factor 
of five, while maintaining 90% efficiency, by setting A ~ 0.1. Selection of a larger 
value of Kq, i.e. deeper wells in the optical lattice, would allow for a larger value 
of A and shorter transfer times. 

4 Conclusions 

We have demonstrated a method for the coherent transfer of ultracold atoms 
from the well of a stationary optical lattice into that of a travelling lattice. The 
effective Hamiltonian for an atom in such a lattice, constructed by adiabatic 
elimination of the internal electronic structure, provides a system of eigenstates 
which determine the center-of-mass dynamics of the atom. We have shown that 
small, harmonic modulations of the lattice amplitude can allow for a STIRAP- 
type subsystem of the Schrodinger equation for this effective Hamiltonian, with 
which transitions between these eigenstates can be induced. 
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time 

Figure 9: The numerical evolution of a state initially localized in the station- 
ary cosine wave of the two-resonanace system (|((/)2(O)|'0(O))P ~ 1), under the 
influence of V{t) for the example described in the text. The probability of a 
transition to any state other than |c) is plotted versus the width of the coupling 
fields (T = (Ji = (72 with evolution under resonant coupling (fii — loi and 
— LO2) plotted with triangles and evolution under near- resonant coupling 
(with Ai and A2 as shown in Figure O plotted with squares. Solid lines are 
the corresponding values predicted by the numerical evolution of the three-state 
model. Figure (b) shows the evolution of the occupation probabilities for the 
case of (T ~ 5000 and resonant coupling (axes on inset Husimi functions, each 
taken at an integer multiple of the period T, are the same as in Fig [7]). 
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A Experimental creation of the effective Hamil- 
tonian 



In this appendix, we show how the effective Hamiltonian in Equation (U) can be 
used to describe an experimental system of lasers impinging on non-interacting 
alkali atoms. This analysis involves the consideration of a two-level subsystem of 
the atom's electronic levels, application of the rotating-wave approximation, and 
adiabatic elimination of the excited level, to obtain a spatially and temporally- 
periodic potential for atoms in the ground state. 

We begin by considering the Hamiltonian of this two level system, in dipole 
interaction with a z-polarized electric field: 



H — Hate 



Hi, 



with 



and 



Hato,n = hujat\e){e\ + ^ (|e)(e| 
Zra 



\9){9\) , 



(55) 

(56) 

(57) 

where hujat is the energy spacing of the two levels, Px is the atomic momentum 
operator in the x-direction, and d = (e|c?z|g) = {g\dz\e) is the dipole matrix 
element coupling the ground state \g) to the excited state |e). The total electric 
field amplitude Ez{x,t) is assumed to be the superposition of the electric fields 
due to N lasers, all polarized along the z direction, so that 



H,nt = -dEz{x,t){\e){g\ + \g){e\) , 



N 



Ez{x,t) = ^£;(^')cos 



(fcz 



Sk.j . , Sujj . 



(58) 



where i?^^-* is the amplitude of the jth laser, is a positive reference frequency 
and fci its corresponding wavevector, aj can be ±1, and 6kj = Sujj/c (the 
usefulness of this form will be evident below) . We can then write 



E,{x, t) = A{x, t) e"''^^' -I- A*{x, t) e 



with 



t)^2^ -IT- exp 



(fc 



Ski , Suji 



(59) 



(60) 



Under a time-dependent unitary transformation of the Schrodinger equation, 
the Hamiltonian transforms like 



H 



BU 

UHU^ + ih—U'' . 
at 



(61) 



Using the unitary matrix 



U = exp |e)(e| t] 



(62) 
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to transform to the rotating frame of the laser leaves the Hamiltonian as 

H = hA\e){e\ + ^i\e){e\ + \g){g\) 

-dE,{x,t){\e){g\e'^-' + \g){e\e-^^-') , (63) 

where A = LUat — i^l is the detuning of the reference laser frequency from the 
atomic transition. 

Let us now make the rotating wave approximation by inserting the form of Ez 
in Eq. (|59p and neglecting terms with high-frequency exponential dependence 
(i.e. e^*^'^^*). The Hamiltonian then takes the form 

H ^ hA\e){e\ + ^(|e)(e| + \g){g\) ^ diAix,t)\e){g\ + A^x,t)\g){e\) . (64) 

Writing an arbitrary state \^) — 'ipg(x,t)\g) + 'ipeix,t)\e) , the Schrodinger equa- 
tion can be written 



^r^if = -dAix,m,+ [nA--—'li... (66) 

Adiabatic elimination of the excited state is performed by assuming that the 
detuning of the laser A is large enough to allow us to neglect the time and space 
derivatives of the excited state. Thus, atoms prepared in the ground state will 
remain there and we are left with an effective Hamiltonian for their evolution: 

'^^-^^ff^s ' HA ■ 

The particular form of A{x,t) will depend on the choice of lasers. A pair 
of counter-propagating lasers with equal carrier frequencies (iJ^^^ — E^"^^ = E\ 
6uji = 5ljJ2 = 0; iTi = —a2\ (/"i = 02 = 0) will produce a time- independent, 
periodic potential, i.e 

Astandix, t) = I (e''^^" + e-"^--^) ^ \A,tand{x, t)\^ ^ ^ C0s{2kLx) , 

(68) 

where we have neglected constant terms. Similarly, two counter-propagating 
lasers with slightly offset frequencies (E^^^^ — E'-^'> = E ; 5uji ~ —8102 = 5u) ; 
(7i = —(72; 4>i = 4>2 = 0) will produce a travelling periodic potential: 

\Atrav{x, t)\^ ^ — COs{2kLX ~ Sujt) . (69) 
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If we combine these two pairs of lasers, we create an effective potential with 
the desired terms of Eq. (|4]) , namely 

H= — — [cos{2kLx) + cos(2fcLa; - 5uj t)] . (70) 

It is clear, however, that the |^(a;,i)p for such a system will also contain un- 
wanted cross-terms, which we have neglected in writing Eq. (j70p . In order to 
minimize the effect of these cross terms, we offset the carrier frequency of the sec- 
ond pair by some amount At^ (£^(1) = £'(2) = £(3) = £(4) = E ; &ujx = &uji = ^ 
; (5w3 — Ao; + iuj ; = Aui — Slu ■ ai = = —02 = —174; and 0i — Vi), 
where ojl » Alj >> 5lo. This yields, 

^2 

I Au;o-res(a;, P ~ — [cos(2fcLx) + cos(2fcia; - Sut) 

+ cos(Afca;) cos(Aajt) + cos(2fcia;) cos(Awi)] , (71) 

where Afc = l^oj/c and we have retained only the highest-order terms in the 
frequencies and wavevectors (e.g. Slo is neglected in the presence of Aw). The 
last two terms in this equation present high-frequency oscillations, depending 
on the particular value of IS.L0. As a concrete example, we can consider a system 
of cesium atoms. In references (22, 23], the laser light was detuned by A 
10"'^"'^ Hz from the D2 line {ojl ^ lO^^Hz) and a modulation of 5lo ~ lO^Hz was 
applied to the standing lattice to affect travelling terms in the effective potential. 
Therefore, an offset of the carrier frequency for the second pair of lasers in the 
hundreds of MHz will satisfy A >> Acj » 6uj, and allow one to safely neglect 
the last two terms in the square brackets of Eq. (TfT]) 24] . 

In order to obtain the Hamiltonian in Eq. ((4]), we change to dimensionless 
variables {p' , x\ H' , t' , uj') as follows. Let p' — , x' = 2k]^x, H' = j^H, 

t' = 8u!rt, and oj' = -^^Suj where the recoil frequency of an atom is ujr — -j;^- 
The Hamiltonian in Eq. ([70)1 then takes the form 

H' = {p' f + K [cos {x') + cos (x' + uj't)] (72) 

where k = — ■ Removing the primes, we obtain the desired Hamiltonian. 

It is important to note that in these dimensionless units, changes in momentum 
due to the interaction of an atom with the lasers are integer- valued. Moreover, 
experimental techniques allow for the preparation of atoms in a very narrow 
range of momentum values about zero [23, l23|. Therefore, the eigenvalues of 
the momentum operator will take only integer values, i.e. p\n) = n\n) with 
n e Z. 

B "Evolution" of a Floquet Hamiltonian 

In this appendix, the {t, t') formalism due to Peskin and Moiseyev [2^, [2^ is 
used to justify the time-parametrization of ki and K2 in the model Hamiltonians 
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in Eqs. ([^ and ([?T|) . These models are each subsystems of a Floquet Hamilto- 
nian which was constructed under the assumptions that ki and K2 were constant 
and the Schrodinger equation was time-periodic. The subsequent parametriza- 
tion of such a system by non-periodic functions of time therefore requires a 
more rigorous explanation. Here, we show that a physical system represented 
by a Hamiltonian with both periodic and arbitrary dependence on time, can be 
associated to Floquet-like Hamiltonian in an extended Hilbert space where the 
periodic time-dependence has been reduced to dependence on a coordinate. This 
Hamiltonian is termed "Floquet-like" because its dependence on the coordinate 
time is identical to a Floquet Hamiltonian's dependence on time. The remain- 
ing arbitrary time-dependence of the Floquet-like Hamiltonian determines, via 
the Schrodinger equation, a dynamics in the extended space from which the 
dynamics of the original system can be recovered. 

Consider the Schrodinger equation for a time-dependent Hamiltonian 



d 



(73) 



where x can be considered a single spatial coordinate or a set of coordinates 
and h has been set to unity by non-dimensionalization of the variables. We will 
associate to H{x\t) a Hamiltonian of one more coordinate Hp{x^t' ]t) which is 
a Hermitian operator in a larger Hilbert space, extended to include this new 
coordinate t' . The relationship between the two Hamiltonians is defined by 



HF{x,t';t) = H{x,t';t) - i 



d_ 

dt' 



(74) 



with 



H{x,t'-t)\t'=t^ H{x-t). 



(75) 



Clearly, H{x\ t) does not uniquely determine Hf{x, t'; t). The time-evolution of 
a state 'tjj{x, t'\ t) in the extended space is governed by the Schrodinger equation 



i— V'(a;, i'; t) = Hf{x, t'; t)'ip{x, t'; t) , 



(76) 



which can also be written 



= H{x,t'\t)-ip{x,t'\t) . 



(77) 



If we take this equation at the cut t' = t, it becomes 



[{l,{x,t'-t)\t,=t\ = H{x;t) [i,{x,t'-t)\t'^t] 



(78) 



where we have used the identity 



^ [i^{x,t'-t)\t>=t] = 



(79) 
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Comparing Eqs. ((75|) and (I75|) . we see that the evolution of a state in the 
original system can be determined by evolution in the extended system using 

lb{x;t)^^{x,t';t)\t'=t, (80) 

and provided the same initial condition 

^(x,t';t)\t'=t=o = ^{x,0) (81) 

is used in each space. 

We now apply this formalism to STIRAP transitions in the two-resonance 
Hamiltonian. The evolution plotted in Figure [SJd was performed by numerical 
integration of the Schrodinger equation, using the Hamiltonian 

H{x; t) — ~'q~^ + '*o [cos a; + cos(a:: — UJot)] 

+ Xcosx[Ki{t) cos(r2it) + K2{t) cos{n2t)] , (82) 

where ki and K2 were given Gaussian time-depenence in order to affect the 
STIRAP-like transition. The solid lines plotted in Figure [9^, were determined 
by evolution of a Schrodinger equation using the time-parametrized three-level 
model in Eq. (j51|) . Using the above analysis we can show that, modulo the 
perturbation theory approximations, these two methods of time-evolution are 
equivalent. We define a Hamiltonian in the extended space 

H{x,t'; t) = + '*o [cosx + cos(x - loqI')] 

+ Xcosx[Ki{t)cos{nit') + K2{t)cos{n2t')] , (83) 

which satisfies Eq. (|75p for the two-resonance Hamiltonian and has the prop- 
erty that functions periodic in time are now functions of the extra coordinate, 
while the amplitudes of the modulations are functions of the usual time param- 
eter. The full Hamiltonian in the extended space HF{x,t':t), defined by Eq. 
((74|) . has the same dependence on t' that the Floquet Hamiltonian in Eq. (I29p 
has on time t. Therefore, the entire perturbation analysis performed on the Flo- 
quet Hamiltonian in Section [3] would proceed in identical fashion on Hp{x, t'; t), 
yielding a time-parametrized three-level model. If Eqs. ((80)) and ([81]) are sat- 
isfied, the "time-parametrized" Floquet Hamiltonian can be used to determine 
the physical evolution. 
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